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Automatic,  Continuous  River  Stage  Measurement 
with  a  Millimeter-Wave  FM-CW  Radar 

NORBERT  E.  YANKIELUN  AND  MICHAEL  G.  FERRICK 


I^ml0Duc^0N 

Radar  sjrstems  are  cximmonly  used  in  hydrologi¬ 
cal  and  other  geophysical  applications,  particularly 
for  measuring  ice  thickness  on  freshwater  rivers  and 
lakes.  The  majority  of  the  effort,  to  date,  has  been 
with  radars  operating  at  frequencies  from  less  dian  1 
to  12  GHz  that  are  able  to  resolve  ice  thicknesses  in 
the  range  of  10  to  20  cm  (Page  and  Ramseier  1975, 
Venier  and  Cross  1975,  Chudobiak  et  al.  1978,  Batson 
et  al.  1984,  Aroone  1991).  Additionally,  radars  have 
been  used  for  other  geophysical  remote  sertsing  ap¬ 
plications  (Wittmann  and  Stoltenberg  1981)  and  for 
snow  stratification  investigations  (EUerbruch  and 
Boyne  1984,  Gubler  and  Hiller  1984).  iladars  are  also 
commonly  used  as  altimeters,  capable  of  measuring 
the  range  from  an  airaaft  to  tlw  ground  (Skohuk 
1980). 

Recently,  an  airborne  millimeter-wave  (MMW) 
FM-CW  (frequency  modulated-continuous  wave) 
radar  system  has  b^  developed  that  can  do  high- 
resolution,  continuous  profrling  of  river  and  lake  ke 
(Yankielun  1992,  nkielun  et  al.  1992, 1993).  It  can 
resolve  distance  in  air  to  as  little  as  1.1  cm,  and  can 
resolve  freshwater  ice  sheet  thicknesses  to  a  mini¬ 
mum  of  3  cm  ±10%. 

Accurate  river  stage  measurement  is  critical  to 
flood  hazard  and  water  resource  assessment.  In  par¬ 
ticular,  the  importance  of  river  surging  on  the  pro¬ 
cesses  of  dynamic  ice  breakup  and  ice-iiuluced 
flooding  is  widely  recognized.  However;  stage  mea¬ 
surement  on  rivers  that  experience  dyiramic  break¬ 
up  can  be  difficult.  Ice  rubble  piles  on  the  batrks, 
darkr\ess  and  ice  scour  of  the  bed  frequently  stop 
measurements  by  survey  teams  and  instruments  re¬ 
quiring  that  cables  be  placed  m  the  river.  It  is  clear 


that  measurement  systems  ttiat  rue  ireadily  installed 
and  not  affected  by  light  and  ice  motion  conditions 
are  needed.  River  stage  measurement  with  high- 
resolution  MMW  FM-CW  radar  provides  a  remote¬ 
ly  mounted  and  monitored,  continuously  record¬ 
ing  alternative  to  the  established  techniques.  Radar 
measurements  from  above  the  river  keep  personnel 
safe  from  the  physical  hazards  of  dir^  observa¬ 
tion  and  eliminate  data  reliability  problems  caused 
by  scour.  As  an  economy,  the  radar  system  used  for 
stage  measurement  is  al^  useful  for  measuring  ice 
thickness  (Yankielun  et  al.  1993). 


FM-CW  RADAR 


In  an  FM-CW  system  (Fig.  1),  the  output  of  an 
MMW  linear  Voltage  Controlled  Oscillator  (VCO)  is 
transmitted  toward  ffie  target.  The  energy  reflect 
from  the  target,  delayed  by  the  round-trip  propaga¬ 
tion  time  Tip,  is  mixed  with  a  sample  of  the  VCO 
output.  The  difference  frequency  F,  is  proportional 
to  the  distarKe  to  the  target  and  can  be  determined 
using  spectral  analysis  techiuques. 

Distance  d  in  meters  is  calculated  in  temts  of  fre¬ 
quency  according  to  the  relation 


d  -  (*»wp)^ 

"2(BW)(n.i,) 


(1) 


where  F,  =  difference  frequency  due  to  radar  re¬ 
flection  (Hz) 

tntp  -  FM-CW  sweep  time  (s) 

c  =  velocity  of  light  in  vacuum  (m/s) 

BW  =  FM-CW  swept  bandwidth  (Hz) 

^•ir  -  index  of  refraction  of  air  =  1. 


(p.'oportlond  to  mdtr  rtng*)  Receivtd  Signal 


Figure  1.  FM-CW  radar  si/stem. 


Figure  2  shows  a  processed  PM-CW  radar  return 
from  a  water  surface^  where 

j  (11.4  kHz) (0.066s) (3xl0®m/fa) 
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The  broader  the  swept  bandwidth  is,  the  greater 
is  the  distance  resolution.  For  the  system  discussed 
here,  the  bandwidth,  as  available  horn  the  sweep  os¬ 
cillator,  is  26.5  to  40  GHz.  There  is  a  tradeoff  between 


where  Frl  B  difference  frequency  due  to  air/ice  in¬ 
terface  reflection 

Fft^  =  difference  frequency  due  to  ice/water 
interface  reflection 

«i£e  =  index  of  refraction  of  freshwater  ice  = 
1.78  (Gumming  1952). 

The  significant  contrast  between  indices  of  re¬ 
fraction  of  air  and  water  causes  radar  scans  obtained 


range  and  resolution  as  there  is  a  limit  on  the 
number  of  time  series  samples  that  can  be  taken 
during  a  given  sweep  at  a  set  sample  rate.  The 
greater  the  maximum  radar  range,  the  fewer 
samples  that  can  be  allocated  to  a  small  he- 
quency  segment  of  interest.  In  this  application, 
the  20^  digitized  time  series  samples  of  each 
scan  are  transformed  to  a  power  spectrum  con¬ 
sisting  of  1024  discrete  frequencies  from  0  to  the 
one-half  the  sample  rate;  the  range  resolution  is 
approximately  1.1  cm. 

When  a  sheet  of  cold,  dry  heshwater  ice  is 
encountered  by  the  radar,  two  pulses  are  re¬ 
flected,  one  eadi  from  the  air/ ice  and  ice/water 
interfaces.  The  ice  thickness  t|  in  meters  is 
found  from  the  separation  of  the  two  difference 
frequencies  according  to  the  relation 
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Figure  2.  Typical  frcqueacy-transfomed  FM~CW  radar  sig- 
(2)  nal  reflected  from  an  air/water  interface  with  radar  range  ««* 
dicated. 
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from  warm  and  wet  ice,  saturated  snow  cover,  sur¬ 
face  melt  pools  or  open  water  to  exhibit  a  single  re¬ 
flection  pulse,  since  virtually  all  of  the  transmitted 
energy  is  reflected  at  this  interface.  The  minute  com¬ 
ponent  of  transmitted  energy  coupled  through  this 
interface  is  predominantly  absorb^  by  the  wet  me¬ 
dium. 


SYSTEM  DESCRIPTION 

The  system  used  for  this  experinaent  (Fig.  3)  is  a 
modihed  version  of  a  successfully  deployed  MMW 
FM-CW  radar  system  (Yankielun  1992,  Yankielun  et 
al.  1992).  Here,  the  system  consists  of  an  HP  3314A 
function  generator  programmed  to  produce  a  0-  to 
10-V  linear  ramp  with  a  period  of  0.056  seconds  that 
is  used  to  drive  the  MMW  VCO.  A  second  output  on 
the  signal  generator  provides  a  pulse,  synchronized 
with  the  start  and  end  points  of  the  linear  ramp,  to 
trigger  an  HP  3660A  Dynamic  Signal  Analyzer  and 
the  12-bit  analog-to-digital  converter  (A/D)  internal 
to  the  laptop  computer.  Data  acquisition  software 
controlling  the  A/D  converter  is  programmed  to 
count  the  synchronization  pulses  from  the  HP 
3314A  and  digitize  every  radar  scan  into  204S 
time  series  samples.  Here,  n  is  chosen  to  permit  ac¬ 
quisition  of  one  scan  every  1, 10  or  60  seconds,  as 


HP3660A 
Dynamic  Signal 

Laptop  Computer  Analyzer 


Figure  3.  Components  of  the  MMW  FM-CW  radar 
system  used  here. 


appropriate  to  fully  characterize  an  event.  Data  are 
sampled  at  31.25  kHz,  providing  a  maximum  radar 
range  of  11.46  m  in  air.  ^ch  digitized  scan  is  sequen¬ 
tially  stored  in  binary  format  on  the  laptop  com¬ 
puter  hard  disk  for  later  Ehgital  Signal  Processing 
(DSP)  and  analysis.  A  time  stamp  and  index  number 
ars  stored  for  each  acquired  scan  in  ASCH  format  in 
a  second  data  file.  The  time  stamp  data  are  used  in 
conjunction  with  the  radar  scan  data  to  precisely 
time  events.  In  parallel  with  this  data  acquisition 
operation,  the  HP3660A  Dynamic  Signal  Analyzer 
is  used  to  monitor  the  process  and  display  real-time 
ranging  information. 

The  radar  front-end — consisting  of  the  VCO, 
waveguide  components,  transmit  and  receive  an- 


Figure  4.  MMW  FM-CW  radar  front-end  suspended 
by  a  cable  from  a  window  in  the  covered  bridge  over  the 
Connecticut  Riixr  at  Windsor,  Vermont.  Initial  range 
from  the  radar  to  the  river  surface  was  approximately  10  m. 
The  slightly  ojf-nadir  antenna  look  angle,  while  not  opti¬ 
mum,  did  not  appreciably  affect  radar  data  acquisition. 
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Relative  Power  (dB)  Scan  Number 


Fi;;;urc  5.  Radar  syatctii  control  and  data  acquisition  components  in  operation  on  the  covered  brid<fe  over  the 
Connecticut  River  at  Windsor,  Vermont.  Vi-dhie  arc  the  laptop  computer  (upper  left),  dynamic  signal  analyzer 
and  sjocep  ;,;eneralor  (lower  left)  and  a  haltery  box  power  supply  for  the  radar  front-end  (lower  ri^hl). 


Frequency (kHz) 


Figure  b.  Relationship  of  a  monochro¬ 
matic  multi  scan  spcctroc^raphic  dis¬ 
play  toil  sin^^ie-scan  poawr  spectrum 
of  radar  data. 
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tennas,  mixer  and  audio  amplifier — is  mounted 
with  the  antennas  pointing  at  a  normal  incidence 
from  a  fixed  point  overlool^g  the  river  (Fig.  4).  The 
remaining  control  and  acquisition  components  are 
remotely  located  and  connected  to  the  front-end  via 
power  and  signal  cables  (Fig.  5).  The  system  is  pow¬ 
ered  by  a  portable,  gasoline-driven  electric  genera¬ 
tor. 

After  a  survey  is  completed,  the  raw  time  series 
radar  data  are  processed  and  displayed  using  a 
Macintosh  II  computer  equipped  with  a  Spectral  In¬ 
novations,  Inc.,  DSP  coprocessor.  Each  radar  scan 
was  digitized  to  provide  2048  time  series  samples, 
transformed  into  a  power  spectrum,  processed  with 
a  Harming  window  algorithm  to  suppress  the  efiect 
of  spectral  sidelobes  that  might  o^erwise  mcisk 
lower  level  signals,  and  displayed  in  a  continuous 
spectrographic  form.  In  a  spectrogram,  discrete  sig¬ 
nal  magnitudes  are  represented  by  a  range  of  color 
or  gray  scale.  With  monochromatic  graphics,  this  re¬ 
sults  in  signal  magnitudes  above  a  preset  threshold 
appearing  as  black  and  below  that  threshold  as 
white  (Fig.  6).  The  DSP  software  permits  a  threshold 
level  to  be  set  for  a  clear  differentiation  of  the  radar 
pulse  reflection  from  fhe  river  surface.  Here,  as  a  re¬ 
sult  of  the  selected  threshold  level  intersecting  re¬ 
flection  pulses  at  a  point  somewhat  below  the  pulse 
peak,  the  trace  appears  fairly  broad.  A  mulfi<olor 
spectrographic  display  provides  significantly  great¬ 
er  graphical  resolution  than  is  possible  with  mono¬ 
chromatic  display  by  indicating  intermediate  levels 
of  signal  intensity  with  a  256-shade  color  gradient. 


RESULTS  AND  DISCUSSION 

Experiments  over  open  water  and  an  ice  cover 
were  conducted  to  ascertain  the  viabUity  of  this  tech¬ 
nique  of  remotely  obtaining  river  stage  data. 

Open  water  experiment 

The  radar  system  was  deployed  over  open  water 
from  a  bridge  across  the  Connecticut  River,  near 
West  Lebanon,  New  Hampshire,  to  acquire  river 
stage  data  prior  to  and  during  the  transition  from 
minimum  flow  to  full  turbine  operation  at  the 
Wilder  Dam  hydropower  station,  2  km  upriver.  Un¬ 
fortunately,  the  height  of  this  bridge  above  the  river 
made  it  impossible  to  suspend  the  radar  sensor 
within  its  specified  operational  range  with  the 
cables  available.  However,  a  faint  trace  that  accu¬ 
rately  tracked  the  river  stage  was  visible  in  the  spec¬ 
trogram  of  the  processed  data.  This  trace  is  a  vesti¬ 
gial  artifact  attributable  to  "aliasing"  in  the  digital 


Time  (hr) 

Figure  7.  Spectrogram  of  radar  ranging 
data  indicating  time  and  river  rise  during 
the  experiment  over  open  water.  U.S.  Geo¬ 
logical  Survey  stream  gauging  measurements 
at  the  same  location  are  also  indicated. 


sampling  process  (Oppenheim  and  Schafer  1975, 
Stanley  et  al.  1984).  Aliasing  is  a  well  known  and 
unwanted  manifestation  that  can  occur  when  digi¬ 
tally  sampling  data.  Typically,  it  is  prevented  or  at¬ 
tenuated  to  negligible  levels  by  appropriate  analog 
filtering  prior  to  digital  sampling.  Here,  however, 
aliasing  was  exploited,  in  efiPect,  to  extend  the  range 
of  the  radar  and  provide  an  accurate  representation 
of  river  stage  (Fig.  7).  While  the  spectrogram  trace  is 
not  as  distinct  as  fiiat  obtainable  with  the  seinsor 
placed  within  the  specified  operational  rango,  it 
closely  tracks  the  data  from  die  gauging  station. 
TVace  clarity  can  be  sigriificantly  improved  when  the 
system  is  operated  within  its  specified  range  limits. 

Ice  cover  experiment 

The  radar  system  was  deployed  from  a  bridge 
over  the  Connecticut  River  on  6  March  1992,  at 
Windsor,  Vermont,  to  capture  the  shape  of  a  large 
water  surge  released  from  the  Wilder  Dam  upriver 
during  a  study  of  controlled  ice  breakup  (Ferrick 
and  Mulherin  1989).  Here,  the  radar  sensor  was 
placed  within  its  operational  range,  permitting  ac¬ 
quisition  of  stage  data  with  improved  trace  clarity. 
A  spectrogram  of  the  raw  stage  data  is  shown  in  Fig¬ 
ure  8a.  Simultaneously  with  the  radar,  a  survey  crew 
acquired  stage  data  at  the  same  location  with  a  level 
and  stadia  rod. 

The  power  spectrum  data  for  each  radar  scan 
were  furtlier  processed  to  precisely  derive  the  radar 
range  from  the  location  of  the  pulse  with  the  maxi- 
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Figure  8.  Spectrograms  of  radar  ranging  data  indicating  time  and  river  rise  during  the  experiment  over  an  ice  cover. 


mum  magnitude  in  each  scan.  The  location  of  the 
pulse  wi^  the  greatest  magnitude  in  each  radar 
scan  represents  the  range  to  the  target  that  provides 
the  greatest  reflectivity.  Here,  there  are  three  pos¬ 
sible  targets  to  which  the  maximum  magnitude 
pulse  In  a  scan  may  be  attributable;  an  air/ice  inter¬ 


face,  an  ice/water  interface  or  an  air/water  inter¬ 
face.  It  has  been  shown  (Yankielun  1992)  that  in  the 
case  of  cold,  dry,  freshwater  ice,  there  may  be  as 
much  as  a  ±15  dB  difference  between  the  magnitude 
of  the  air/ice  and  ice/water  reflected  pulses  within 
a  given  scan,  mainly  attributable  to  variations  in  the 
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surface  and  volume  scattering  of  the  radar  sig¬ 
nal.  Therefore/  in  any  given  scan  taken  over 
cold,  dry  ice,  the  air/ice  or  ice/water  pulse 
may  be  identified  by  die  peak  finder  algorithm 
as  the  pulse  with  the  maximum  magnitude. 

This  results  in  a  possible  scan-to-scan  variation 
of  range  proportional  to  the  thickness  of  the  ice 
sheet.  Figure  8b  shows  the  result  of  processing 
with  the  maximum  peak  locating  algorithm. 

Data  from  the  survey  measurements  are  over¬ 
laid  for  comparison. 

To  smooth  the  data,  a  miming  10-point 
averaging  algorithm  was  applied  to  the  maxi¬ 
mum  peak  magnitude  data,  wifii  the  results 
shown  in  Figure  8c.  Here,  the  overlaid  survey 
data  aoc  uately  track  the  radar  data  after  1300 
hours.  Prior  to  this,  the  averaged  radar  data  are 
either  systematically  higher  of  shift  irregularly 
about  the  survey  data.  The  difference  between 
the  radar  and  the  survey  measurements  is  caused  by 
shifting  ice  plates  in  the  near  bank  area  as  they  inter¬ 
act  with  each  other  and  the  river  bed.  Loca^g  the 
radar  system  further  from  the  bank  over  smooth 
and  free  floating  sheet  ice  would  provide  accurate 
river  stage  measurement. 

At  0730  hours  on  30  March  1993,  the  radar  was 
again  deployed  over  the  Cormecticut  River,  at  Wmd- 
sor,  Vermont,  to  obtain  stage  data  during  a  natural 
breakup.  The  ice  had  broken  into  mbble  and  was  al¬ 
ready  moving  when  the  radar  system  was  set  up 
and  activated,  The  stage  data  were  acquired  every 
20  seconds;  Figure  9  shows  the  resulting  stage  radar 
spectrograph.  The  wide  variation  seen  in  the  left  half 
of  the  spectrograph  trace  results  fiom  the  variation 
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Figure  9.  Stage  data  obtained  during  30  March  19S2 
showing  the  transition  from  a  moving  rubble  ice  cox>er 
to  open  water  on  the  Connecticut  River  near  Windsor, 
Vermont. 


Figure  10.  Comparison  of  radar  and  surveyed  river  stage  data 
obtained  on  the  Connecticut  River  near  Ascutney,  Vermont, 
on  30  March  1992. 


in  rubble  height  above  the  water  surface,  and  indi¬ 
cates  t!  le  trend  in  river  stage.  During  the  second  half 
of  the  event,  the  ice  mbble  had  passed  from  under 
the  radar  and  the  resulting  scan-to-scan  consistency 
is  indicative  of  a  response  fiom  open  water.  No  sur¬ 
veyed  stage  data  were  available  for  comparison  to 
the  radar  measurements  during  this  portion  of  the 
breakup  event  owing  to  the  hazmxi  presented  by  the 
moving  ice.  The  stage  increase  following  the  pas¬ 
sage  of  the  ice  continued  over  the  next  hour,  but  at  a 
diminishing  rate.  During  this  same  period,  the  stage 
was  steady  at  a  stream  gauge  23  km  upstream  in 
open  water,  and  at  a  survey  station  8  km  down¬ 
stream  where  the  river  had  a  stationary  ice  sheet. 
Together,  these  data  indicate  that  the  rising  stage  at 
Wmdsor  was  local,  resulting  from  the  ice  stopping 
and  a  backwater  developing  behind  the  newly 
formed  jam. 

At  0930  hours  on  30  March  1993,  the  radar  sys¬ 
tem  was  redeployed  from  another  bridge  over  the 
Connecticut  River  at  Ascutney,  Vermont,  approxi¬ 
mately  8  km  downstream  of  the  Windsor  site.  Figure 
1C  compares  the  radar  spectrogram  with  surveyed 
stage  data  at  this  site.  The  agreement  between  these 
data  records  is  excellent  for  ^e  entire  period  of  mea¬ 
surement.  The  river  stage  increased  gradually  for 
several  hours  and  then  increased  abiuptly.  The 
rapid  rise  beginning  at  1358  was  caused  by  the  ar¬ 
rival  of  a  flow  surge  initiated  by  sudden  ice  motion 
upstream.  The  local  ice  motion  began  at  1407  and 
stopped  at  1441,  just  prior  to  the  river  stage  peak. 
After  the  peak,  the  stage  rapidly  diitdnished,  reach¬ 
ing  a  new  steady  state  prior  to  1600.  'fhe  shape  of 
this  portion  of  the  stage  graph  is  cha-acteiistic  of  a 
"drainout,"  indicating  that  the  source  of  the  surge, 
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Icxral  channel  storage  because  of  the  ice,  was  hnite 
and  had  been  depleted. 


CONCLUSIONS 

Few  river  stage  measurement  methods  are  avail¬ 
able  that  can  be  rapidly  deployed  and  are  reliable 
during  cdl  conditions,  and  moving  river  ice  poses  a 
great  hazard  to  those.  These  experiments  over  both 
open  water  and  an  ice  cover  demonstrate  that  an 
automatic,  stand-alone  MMW  FM-CW  radar  sys¬ 
tem  can  be  rapidly  installed  to  provide  continuous, 
around-the-clock  stage  data  of  accuracy  comparable 
to  those  acquired  by  a  survey  team  lifted  to  day¬ 
light  operation.  Since  the  system  can  be  remotely 
mounted  and  monitored,  persotmel  are  not  jeopar¬ 
dized  by  the  hazards  involved  with  direct  measure¬ 
ment.  ^dar  stage  measurement  from  above  the 
river  eliminates  the  data  reliability  problems  of  sys¬ 
tems  that  require  the  placement  of  cables  in  the  river, 
which  are  subject  to  scour.  The  radar  system  is  also 
economical  because  it  is  identical  to  that  used  for  ice 
thickness  measurement. 

Our  initial  experiments  show  promise  and  have 
pointed  out  several  modifications  for  future  mea¬ 
surements.  A  mounting  bracket  will  be  installed  to 
eliminate  any  wind-induced  movement  of  the  radar 
front-end  over  the  river  surface.  A  better  maintained 
nadir  look-angle  will  ensure  that  the  maximum 
amount  of  transmitted  power  will  be  reflected  di¬ 
rectly  back  to  the  radar,  thus  improving  the  signal- 
to-noise  ratio  of  the  system.  Additionally,  position¬ 
ing  the  radar  more  centrally  in  the  river  over  free- 
floating  sheet  ice  rather  than  closer  to  the  bank 
avoids  shore  and  grounded  ice  sheet  effects  and  im¬ 
proves  the  accuracy  of  river  stage  measurements. 
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